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Abstract 

We show that the reservoir influence can be modeled by an infinite array of beam splitters. 
The superposition of the input fields in the beam splitter is discussed with the convolution 
laws for their quasiprobabilities. We derive the Fokker-Planck equation for the cavity field 
coupled with a phase-sensitive reservoir using the convolution law. We also analyse the 
amplification in the phase-sensitive reservoir with use of the modified beam splitter model. 
We show the similarities and differences between the dissipation and amplification models. 
We show that a super-Possonian input field cannot become sub-Poissonian by the phase- 
sensitive amplification. 


1 Introduction 

A cavity with imperfect mirrors not only lets the cavity field out but also allows the field outside 
to leak into the cavity so that the reservoir surrounding the cavity gradually influences the cavity 
field. The mirrors in the cavity can be considered as beam splitters. An arrayx consisting of an 
infinite number of beam splitters can model a reservoir coupled to the cavity field [1]. 

It is convenient to utilise quasiprobability distributions such as the Q function to describe 
states of quantum-mechanical systems in phase space. We briefly show that the fields at the 
output ports of the beam splitter can be expressed by the convolution of the quasiprobabilities of 
the input fields. We then use the convolution relation to derive the Fokker-Planck equation for 
attenuation of the cavity field coupled with the phase-sensitive reservoir. 

It is known in quantum mechanics that an amplification process is inevitably accompanied 
by the increase of the quantum noise in the system. In other words the amplification degrades 
an optical signal and rapidly destroys quantum features that may have been associated with the 
signal. The nature of the amplifier affects the physical properties of the amplified states of light. 
The phase-sensitive amplifier is conceptually based on the establishment of squeezed light and 
enables a squeezed input to keep the property for a gain larger than the cloning limit [2]. 

The amplification process can also be modeled by an array of beam splitters, where the beam 
splitters are somewhat modified from the usual sense. We find the convolution relation for this 
modified beam splitters and the Fokker-Planck equation is derived from it. We then formaly solve 
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the Fokker-Planck equation for the phase-sensitive amplifier for an arbitrary input field and study 
the photon statistics of the amplified field. 

2 Phase-Sensitive Attenuation 

Consider that two fields at the two input ports of a lossless beam splitter are superposed. For 
convenience we call one input field the signal and the other the noise (Fig. 1). The input signal 
mode b with its annihilation operator b is superposed on the noise mode a with its annihilation 
operator a by the beam splitter whose amplitude reflectivity^ is r = sin 9 and transmittivity 
t = cos 9. The two output field annihilation operators c and d are related to the beam splitter 
input fields by the transformation using the beam splitter operator B [3]. 



A 

a 


Fig. 1 Beam splitter with the signal input in mode b and the noise field input in mode a. 

A field state can be represented in phase space by quasiprobabilities. Let us choose to use the 
positive P-representation to describe the fields. The density operators p a and p b for the noise and 
signal fields are then written with the positive P-representationP a>b (a, 7 ) as 

A = /<farf , 7 4^ |o)( 7 *l. i = a,b ( 2 ) 

The fields p c and p d at the output ports are calculated by the beam splitter transformation: 
•BpaPb — * Pc Pd- With use of the transformation matrix Eq.(l) we find that the positive P- 
representation for the field at the port d is found in the form of the convolution relation: 

Pd(M) = ^J^ a d2 i f r ~ . ~ > ( 3 ) 

The Q function is another quasi probability function and is well-defined even for the nonclassical 
state. The positive P-representation is defined in four-dimensional space. The Q function, which 
is defined in two-dimensional space, is therefore sometimes easier to treat, so we will extend the 
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convolution law for the uses of the Q function. The positive ^-representation may be defined as 
the Fourier transform of the characteristic function. The characteristic function C (p) (£) is related 
to the characteristic function C^(0 for the Q function as C (p) (0 = C (q) (f) exp(|f| 2 ) . By using 
the convolution theorem, we can factorise the inverse Fourier transform of the convolution law (3) 
as 

Ci p) ( f) = C<-Hr()Ci P, (tO . (4) 

Using the relation between the characteristic functions, the convolution relation for the Q function 
is found as . , \ 

QM = ^J d2 « QMQ» • ( 5 ) 

We now derive the Fokker-Planck equation for the phase-sensitive reservoir using the model 
of an infinite array of beam splitters (Fig. 2) [4]. The total duration of time when the field is 
coupled with the lossy channel is denoted by T , the total number of the beam splitters by P, and 
the interval between the adjacent beam splitters by At. The beam splitters are first taken to be 
discrete components, but their number, N = T/Ar is later taken to infinity in order to model a 
continuous attenuating reservoir. Under the assumption that the reflectivity is very small for the 
beam splitter, Eq. (5) is written as 

Q(a) « (1 + R)J (fa QM Q b • ( 6 ) 


N beam splitters 



P out 


Fig.2 The phase-sensitive reservoir modeled as an array containing an infinite number of 
beam splitters. The signal is injected from left and the independent squeezed fields (all with the 
same properties) are injected into the other ports. The transmittivity is considered to be nearly 
unity. 


To calculate the effects of attenuation, we need an expression for the output signal operator in 
terms of the input operators. To simulate an attenuator, we consider the beam splitters forming 
a continuous array by taking the limits P — * oo, At — * 0, and Ft * 0 . These limits cannot 

be taken independently: KAt should be kept constant. Also, the total energy loss within T 
is described by 1 - exp (- kT), where k is the attenuation coefficient, and this loss should be 
equivalent to the beam splitter loss so that R « kAt. 
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Let us define Q(r; a) as the Q function of the signal field incident on the beam splitter at 
time r, Qsq(/9) as the Q function for the noise added to the signal at the beam splitter, and 
Q(t + At; a) as the Q function for the signal leaving from the beam splitter. The squeezed 
thermal fields produced by the independent stationary sources act as noise in our model. 

From Eq. (6), we obtain the relation 


Q(r + At; a) = (1 + R) j d 2 P Q (t; Q s. 


(/?) 


(7) 


where - « a + ^ a — \/R 0. With the Taylor expansion for a function having a complex 

argument: 
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where the real and imaginary parts of a and P are respectively denoted by aj, a 2 and Pi, /? 2 . The 
function Q is the simplified notation of the function Q(t\ a). Substituting Eq. (8) into Eq. (7) we 
obtain 


dQjr] a) 
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Taking the squeezed thermal state as noise we substitute the simple Gaussian integration of 
Qsq{P) into Eq. (9) and obtain the Fokker-Planck equation for the field coupled to a phase-sensitive 
attenuation reservoir: 


dOi^a) 

dr 


5^“' + + 5 (1 + N + + 5 (1 + N _ m) 55|J 


Q{r,a) , (10) 


where N is the mean photon number for the squeezed thermal field and the phase-dependent term 
M is zero when the field is not squeezed [4]. The Fokker-Planck equation is relatively simple and 
observables can be calculated as correlations of the quasiprobability function. 


3 Phase- Sensitive Amplification 


The amplification process can also be modeled by an array of beam splitters similar to that for 
dissipation. In their experiment with phase-sensitive amplification, Ou et al. have a nondegenerate 
parametric amplifier where the signal field is amplified and the idler mode is coupled with the 
squeezed vacuum [5]. As in their experiment a two-mode parametric optical amplifier is modeled 
here by an amplification beam splitter matrix. For a two-mode parametric amplifier the signal 
input b is transformed into the amplified output d with unavoidable noise fit : 
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( 11 ) 
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where g > 1 is the infinitesimal amplification factor. We are going to build a beam-splitter-like 
relation for amplification, and consecutive application of an infinite number of Eq. (11) will give the 
final amplification result. The actual gain G by the amplifier will, thus, be proportional to g . The 
unitary amplification beam splitter operator has been introduced in analogy with the two-mode 
squeezing operator as B\ — exp[i0i(a^6^ + a6)] with cosh$i = yj~g and sinh theta i %/g — "T. 

To analyze the beam splitter transformation for the amplifier, let us assume that the input 
fields are expressed as a weighted sum of diagonal coherent components: 


p in = I d 2 a d 2 p P a (a)P b (/3)|a)aa(a| ® \P) hb {P\ , ( 12 ) 

where P a and l\, are respectively the Glauber P-representations for modes a and b. Tracing the 
output field over mode c, we find the output density operator for mode d: 

p d = J d 2 a d 2 p P.(a)ft(flA,(«)pthI>l(«) ; 6 = isfd^ <** + V9 P (13) 

and Ah is the thermal field density operator for the mean photon number n - g - 1. Even when 
both the signal and the idler fields are in the vacuum state, i. e., a = 0 = 0, the amplifier 
brings noise into the fields. The density operator for the thermal field can be written with its 
quasiprobability Pr(0) as 

Pth = J d 2 (f> Pr(<t>)\4)(4>\ 5 = ^3Y ex P 

By using Eq. (14), we find the Glauber P- representation for the output field as convolution of the 
three P-representations 

PM) = i/dV awft p T (« . (15) 

The inverse Fourier transform of the Glauber P- representation gives the characteristic function 
for the output field in the form of the product of the characteristic functions for the input modes 
a and b and the thermal field. We can simplify this relation using the relation between the 
characteristic functions for the various quasiprobabilities. The Fourier transformation of this 
shows that a modified convolution between Q b for the signal Q function and P a for the noise 
P-representation results in Qd for the output field: 

QM) = - g J d2x a(A)0b ~ " / jf 1 A ) ■ ( 16 > 

The convolution relation for amplification differs from that for attenuation (5) because of the 
unavoidable extra noise due to the thermal field (14). 

Consider an array of N beam splitters which satisfy the transformation relation (11). To 
simulate an amplifier we will take K — ► oo and let the infinitesimal amplification factor for each 
beam splitter be given by g — 1 4* f w 1. After a signal passes through the K beam splitters, it is 
amplified by the factor of G = e lT = (1 + e)*, where 7 is the amplification coefficient. By using 
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the Taylor expansion of Q function (16) to the second-order under the assumption t « 0 (we had 
r « 0 for attenuation) , we obtain the Fokker-Planck equation for amplification: 


dQjr, a) 

dr 
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Q(r,a). (17) 


The Fokker-Planck equation (17) is solved for an arbitrary signal amplified in the phase-sensitive 
reservoir. The Q function corresponding to an arbitrary input field can be written as a weighted 
integral of Gaussian functions: 

Qs{a) = — f d 2 fid 2 uP(n, v) exp[-(au - A) 2 - {a 2 - B) 2 } , (18) 

7 r J 


where A = l(/x + i/), B = j^—v) and P(/x, v) is the positive P-function for the field. It has been 
recently shown that if the initial Q function of the quantum system is (complex) Gaussian, then 
the solution of the Fokker-Planck equation (17) is also Gaussian with time-dependent parameters. 
The Q function (18) is a weighted integral of complex Gaussian functions, so one can obtain the 
time evolution of the input state 


Qamp(o ! > T ) 
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(19) 


where the time-dependencies of the amplification parameters are 

A(t) = Ay/G = Ae^ , Bit) = B\[G , JV,(r) = (N + 1)(G - 1) and M,(t) = M ( G - 1) . 

( 20 ) 

The inverse Fourier transformation of the Q function (19) shows that the characteristic function 
C^JJjp for the Q function of the amplified field is the product of the characteristic function 
for the P-representation of the squeezed thermal field and that for the Q function of the 
amplified signal without noise 


CSt(C) = C<J>(iv'g^7 <)C2">(V5 <) . (21) 


The convolution relation for this relation is then in a form analogous to Eq. (16) for the amplifi- 
cation beam splitter superposition of two input fields. 

The antinormally-ordered moments can be calculated from the characteristic function C^: 


(s n (a') m ) 
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( 22 ) 


We can also consider the phase-sensitive amplifier which can be implemented as a stream of 
three-level atoms in a ladder configuration with equispaced levels injected into the cavity where 
the initial state of the field has been prepared. We denote the population in the uppermost state 
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by Paa, the population in the lowest state by p cc and the coherences between them by and p, 
The atomic coherences p & c and pea bring about the phase-sensitive effect in the two-photon linear 
amplifier. The parameters N and M can then be represented by the atomic variables 

Ni(r) = — — — (C - 1), = -J^L_(C - 1) . (23) 

Paa — Pcc Paa Pcc 

The normally-ordered photon number variance, : (An) 2 where (An) 2 = (n 2 ) — (n) 2 , measures 
the deviation of the photon number fluctuations from the Poissonian photon statistics. The Pois- 
sonian field has the normally-ordered photon number variance zero, while the quantum mechanical 
sub-Poissonian field has it less than zero. It is larger than zero for the noisy super-Poissonian field. 
With use of Eq.(22) we find the normally-ordered photon number variance 


: (An ) 2 == G 2 : (An ) 2 : +C, 


(24) 


where the additive noise is 


c = 


(G-l ) 2 

(Paa Pcc) 


(PL + iPad 2 ) + 


G(G ~ ‘W fc'a). - |fel((a 2 ). + <(a') 2 },]. 

Paa Pcc 


(25) 


If the additive noise is negative, the amplified field has less photon number fluctuation than the 
input field. It is clearly seen that if the atomic coherence, p ac , is zero the additive noise is always 
positive. However as the atomic coherence is nonzero we can have the negative noise to enhance 
the signal to noise ratio. 

If each atom injected into the cavity is in atomic coherences we have the relation PaaPcc = iPad 
and the additive noise (25) can be written as 


C = G{G-l) 


c|>> “ 1 + 2 G(G - 1) P “ , fe ; (a t a)i. + (G - 1) ! 

Paa Pcc Paa I Pac I 


Paa 


(Paa Pcc) 


(26) 


c = 2(ata) in -((a 2 ) in + ((a0 2 ) i n), ( 27 ) 

which has to be negative to have the additive noise negative. The bosonic operators a and a f have 
a simple restriction, 2(<Pa) — ((a 2 ) + (( a *) 2 )) ^ — !• ^ thus required that 

— 1 < c < 0 ( 2 8) 


for the noise reduction in the amplified signal. The noise reduction in the photon number fluc- 
tuations seems to be possible if the input field satisfies Eq.(28). However we should not fail to 
notice that the condition (28) is related to the initial photon number fluctuations. Because the 
expectation value of an operator times its hermitian conjugate is again positive, 


c > 


\ (An) 2 , : 

(<da) 


(29) 


It is easily seen from Eqs.(28) and (29) that the input field should be super-Poissouian to have 
a possibility to reduce the photon number fluctuations by the amplification. If the input hvid is 
Poissonian there is no intersection between the two conditions (28) and (29) so that we can say 
that the Poissonian field does not become sub-Poissonian during the amplification. 
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